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Abstract 

We study non-Abelian Chern-Simon BPS-saturated vortices enjoying N = 2 
supersymmetry in d = 2 + 1 dimensions, with generic gauge groups of the form 
U(l) x G", with G' being a simple group, allowing for orientational modes in the 
solutions. We will keep the group as general as possible and utilizing the power- 
ful moduli matrix formalism to provide the moduli spaces of vortices and derive the 
corresponding master equations. Furthermore, we study numerically the vortices ap- 
plying a radial Ansatz to solve the obtained master equations and we find especially 
a splitting of the magnetic fields, when the coupling constants for the trace-part 
and the traceless part of the Chern-Simons term are varied, such that the Abelian 
magnetic field density can become negative near the origin of the vortex while the 
non-Abelian part stays positive, and vice versa. 
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1 Introduction 



Solitons play a crucial role in a vast area of physics ranging from particle physics and cosmology to 
condensed matter physics. Planar physics i.e. in 2+1 dimensions, has radically different properties 
as the spin is not quantized as in 3+1 dimensions giving rise to the anyons among others, objects 
having fractional spin and statistics. This can be realized by the Chern-Simons term which has 
been widely used in e.g. the theory of the fractional quantum Hall effect [T]. Another aspect 
is that the high temperature limit of a four dimensional theory can be described by a three 
dimensional one, where the Chern-Simons term resides naturally. Another interesting feature of 
Chern-Simons theories is that it provides a gauge invariant mechanism of mass generation [2]. 

The most celebrated vortex solution, namely the Abrikosov-Nielsen-Olesen (ANO) vortex was 
found half a century ago [31 H]. This object carries magnetic flux in its interior. Later, similar 
vortex solutions, however in 2+1 dimensions where found with a Chern-Simons term instead of 
a Maxwell term [HIE]. These vortices possess the already mentioned features of fractional spin 
and statistics, viz. they are anyon-like. Furthermore, there exist vortices in both the asymmetric 
phase (like for the ANO vortices) and also in the symmetric phase. The latter do not have a 
topological argument for stability. The vortices with Maxwell or with Chern-Simons terms split 
into three categories depending on the self-coupling of the Higgs field, viz. type I /II vortices or 
the critical BPS saturated vortices [7], where the vortices attract, repel and do not feel any force 
among the selves, respectively. The latter corresponds to some amount of supersymmetry being 
present in the theories at hand. Recently, a fourth type of vortices in the Abelian Chern-Simons 
model has been found, behaving as a type I vortex at small amount of magnetic field and turns 
into type II when the magnetic field piles up repelling further vortices from the clusters [5] . This 
type of vortex was denoted a type III vortex. 

A few years ago, non- Abelian vortices have been discovered 0, HD] , being flux tubes which are 
carrying orientational modes. These models have been extensively studied with the gauge group 
U(N) ~ U(l) x SU(N)/7*n and only recently with generalizations to other groups [TTJ [T2], flB"] . 
In particular the moduli space of these vortices have been studied in detail [5J [TTJ], [J3J HU E2 EH] • 
Good reviews summarizing many results can be found in Refs. [TTJ HU [T9] . 

The first studies of non- Abelian Chern-Simons vortices are made with a simple group, viz. SU (2) 
and SU (N) with fields in the adjoint representation [201 EH [22] and later numerical solutions 
have been found [23]. In Refs. [24], [25] the non- Abelian Chern-Simons vortices have been studied 
with a U (N) gauge group allowing for orientational modes to be present and they identified the 
moduli space of a single vortex solution. Furthermore Refs. [26l [27] [28] have considered packag- 
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ing together the Yang-Mills and the non-Abelian Chern-Simons terms for U(N) gauge groups. 
In Ref. [26] the dynamics of the vortices has been studied and in Ref. [27] in addition to the 
topological charge, conserved Noether charges associated with a [/(l) 7 ^ 1 flavor symmetry of the 
theory due to inclusion of a mass term for the squarks. In Ref. [28] numerical solutions have 
been provided. 

Many related topics can be found in the excellent reviews j29j [30] . 

It is the purpose of this paper to consider a wider class of non-Abelian Chern-Simons vortices 
carrying orientational modes, with the gauge group kept as general as possible, except when we 
will do some concrete numerical calculations. 



2 The model 

Our starting point will be the Yang-Mills-Chern-Simons-Higgs theory. We are considering the 
following M = 2 supersymmetric theory (viz. with 4 supercharges) in d — 2 + 1 dimensions with 
the gauge group G = U(l) x G', where G 1 is a simple group. The bosonic part of the Lagrangian 
density reads 

r^ nr ~„„ — _ ( F a \ 2 — ( F° \ 2 ^LpWP ( A a p) A a — - f ahc A a A b A c I —c^p A f) A 

*-YMCSH - , 2 l/VJ 4e 2 IV J g7r e \ A ^ A P 3/ A p A v A p I g7r e AppvAp 

" T ( Tr ^ " M - 7 ( Tr ^ - " ' (2 ' 1} 

where a = 1, . . . , dim(G"), the index is for the Abelian group and a = 0, 1, ... , dim(G') and we 
use the conventions 

= d„A v - d v A fl + 1 [A» A v \ , V„H = (fy + iA„) H , P M = «9 M + 2 0] . (2.2) 

A^ = A^t a is the gauge potential, is the field strength, is an adjoint scalar field which 
we can take to be real and finally if is a color-flavor matrix of dimension N x Nf of Nf matter 
fields. We will define N = dim (.Re) but for simplicity we choose the representation Rqi as the 
fundamental one of G' . We are using the following normalization of the generators 

*° = , Tr (t a t b ) = -5 ab . (2.3) 

There are four coupling constants entering our game at this point; e G K is the Abelian coupling 
of the Yang-Mills kinetic term (Maxwell), g G K the is the coupling for the semi-simple part 
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of the Yang-Mills kinetic term, which corresponds to G' . k G IR is the Abelian coupling of the 
Chern-Simons term while [i G Z are solely integers to render the non-Abelian Chern-Simons 
action gauge invariant up to large gauge transformations [31]. £ is a Fayet-Iliopoulos parameter. 
Finally, m is a mass matrix which we will set to zero in this paper. 

The scope of study in this paper will be on the Chern-Simons part of this theory. A detailed 
study of the vortices dependence of the parameters of the model above with also the Yang-Mills 
term in action will be done elsewhere 1 3 2 1 . 



3 Non— Abelian Chern-Simons— Higgs theory 



Now let us take the limit e - 
scalar field <fi: 

Air 



oo, g — > oo, m = and and in turn integrate out the adjoint 



Tr (HHH a ) , 



Air 1 



This leaves us with the non-Abelian Chern-Simons theory 



[Tr(Mt) _£] 



(3-1] 



C 



CSH 



.JL^p I A a a d v A a D - -f abc A a u A b A c 

ct. \ p v p t^J pup 



1 

3* 



|-e^ {Ald v A° p ) + Tr (V^ (V»H) 



- 4vr 2 Tr 



— (Tr (HH*) -£) + -Tr (HHh a ) t a \ H 
Nk v v ' ' fj, v ' J 



(3.2) 



which will be the main focus of this paper. It still enjoys M = 2 supersymmetry and there are 3 
parameters governing the solutions; the Abelian Chern-Simons coupling k and the non-Abelian 
Chern-Simons coupling /i and finally the winding number v = ^ [12]. denotes the greatest 
common divisor (gcd) of the Abelian charges of the holomorphic invariants of G', see [12]. For 
simple groups this coincides with the center as Z no . We will take k > 0. 

There are three different phases of the theory at hand. An unbroken phase with (H) = and 



a broken phase with (H) = y jj- In between there are partially broken phases. We will consider 
only the completely broken phase in this paper. 

The equations of motion are 



pu 



— iTr 



8vr 

— c^F* = -iTi 
8tt ^ 



Hk a V a H - (V a H) ] t a H 
HH°V a H - (V a H)U°H 



(3.3) 
(3.4) 
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Vp^H 



Air 2 
Sir 2 



^ (Tr (Mt) - + ^Tr (m¥) r 



167T 2 



Tr 



H 



' (Tr (Mt) - e) + -Tr (HHh a ) t a 
A 4 



Tr 



^(Tr( M t)_ e)+ ^ Tr(M¥)t , 



HH^ j H 
HHhA t a H . 



(3.5) 



The tension, defined by the integral on the plane over the time-time component of the energy- 
momentum tensor, is given by 



T 



Tr <! \V H 



\T>iH\ + 47r" 



^(Tr(Mt)-e) + ^Tr(ififV)r^| 2 | , 



(3.6) 



which by a standard Bogomol'nyi completion can be rewritten as 

T= J Tr j £> # - z2tt f-^- (Tr - f) + ^Tr (HHH a ) t a ^j H 

F° 2 + zTr / [<9 2 {rfVxH) - d 1 (H^V 2 H)] . 
Jc 



4 \VH\ 



(3.7) 



This leads immediately to the BPS-equations which need to be accompanied by the Gauss law 
being the a = component of the Eqs. (I3.3l) . (l3.4p 



VH = , V H = i2iv (J^- (Tr (M f ) - £) + ^Tr (HHH a ) t a ^j H . 



(3. 



Rewriting the boundary term using the first BPS-equation, we have for the BPS saturated vortices 
the tension 



T 



'2N 



j F° + ^Tr / d 2 (Mt) = 2<l 
Jc z Jc 



(3.9) 



with v being the U(l) winding number. By combining the BPS equations with the Gauss law, 
we obtain the following system 

VH = , (3.10) 

1 f)7T^ 1 f)7T^ 

F° 2 t a = (Tr (Hrf) - Tr (HHH a ) t a + -^Tr (HHH b ) Tr (M* {t\ t b }) t a , (3.11) 



57T 



16tt" 



-Tr (Mt) (Tr ( M t) _ ljV + (Tr ( M tf))' x 



^ - v — , v „ v — , v , ^ v - v _ „ „ ^ . (3.12) 

An interesting comment is that the system only depends on three combinations of the couplings; 
viz. k, 2 , fi 2 and nfi. There are thus only two choices of signs giving different solutions sign(ft) = 
±sign(/i). This system is of a generic character and one can readily apply one's favorite group. 
Setting k = fi, the BPS-equations become 



VH = , V n H 



i2n 

K 



2Tr (HHh a ) t a 



1 

N 



H 



(3.13) 
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which in turn yields the simplified system by combination with the Gauss law 



VH — , F&t 



16tt 2 



Tr (M f {t a , t 13 }) Tr (HHH P ) - ^Tr (HHh a ) 



t a . 



(3.14) 



In the next section, we will consider the cases of G' = SU(N), G' = SO(N) and G' = USp(2M), 
and finally make the corresponding master equations. 



3.1 Master equations 
3.1.1 G' = U{1) x SU{N) 

Considering the case of U{1) x SU(N), the BPS-equations combined with the Gauss law read 



VH = 



rpa +a 



in- 



(Tr (Mt) - e) (Hrf - ^Tr (M^ 



(3.15) 



r 12 L 



87T 

+ — 

/I 2 

87T 2 
N 2 K 2 



N 



HH^ Htf - ^Tr (Mt) 
Tr( M t) (Tr(Mt)-e) 1. + g 



Tr ((Mt) 2 
Tr ((M f ) 



AT2 



(Tr (Mt)) : 



LA? 



In this case, the generic vacuum is given by 



(H) 



-A 



(3.16) 



This vacuum allows for an unbroken global symmetry, the so-called color-flavor symmetry which 
is the global part of the gauge transformation combined with the flavor symmetry. This is of 
crucial importance for having orientational modes in vortex configurations. 

Utilizing the moduli matrix formalism, we can immediately solve the first BPS-equation and 
rewrite the second in terms of the new variables 



H = S^Hoiz) , A a t a = -iS'^dS' , A°t° = -idlogs 



(3.17) 



along with the definitions Q = ouQ', Q' = S'S , ou = ss^ and f2 = H (z)Hl(z). The field-strength 
matrices are 



F? 2 t a = 2S'~ 1 d 



S' , F? 2 t° = -21 N dd log to . 



(3.18) 
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In this U(l) x SU(N) case we can write down the two master equations like 

4tt 2 1 fl 







/ Qn /-1 



NKjJbOJ \u 
4vr 2 1 

+ 



Tr Q Q 



Votl'- 1 - ^Tr | 



yU 2 LU 2 



fW" 1 m fi /_1 - ^Tr (fiofi'^H 



(3.19) 



dd logo; 



4tt 2 1, 
iV 2 K 2 ^ 
4vr 2 1 

+ 



( -Tr ( 



Tr ((Wr- 1 ) 2 ) - 1 (Tr (fi ^ 



(3.20) 



Setting the couplings equal k — //, we can write the U(N) Chern-Simons BPS-equations and 
master equation as simple as 

r 2 



^12* 



57T 



4tt 2 

a [nan- 1 ] = — finfT 1 



/V 2 



(3.21) 
(3.22) 



The boundary conditions for these master equations coincide with the weak coupling solutions 

(ism 



3.1.2 G" = C/(l) x SO{N) and G" = C/(l) x USp{2M) 

Considering now the gauge group G — U(l) x SO(N) and G = C/(l) x USp(2M) on the same 
footing with their corresponding invariant tensor J, which has the properties Jtj = In and 
J T = e J with e = ±1 for SO(N) and USp(2M), respectively. 

The vacuum has the generic form [33J 

(H) = diag(ui,u 2 , • • • ,vn) , vteR+, (3.23) 

however, we will consider the most symmetric vacuum allowing for the global color-flavor sym- 
metry, viz. we will here use (I3.16p . We have the following system which is obtained by combining 
the BPS equations with the Gauss law and applying respective algebras 

VH = , (3.24) 

= ^- (Tr (rat) _ e) (m* - Jt (Mt) T j) + ^ [(rat) 2 _ 7 t ((rat) 2 ) T j , 
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which lead to the master equations 

2vr 2 1 (\ 







Tr [ Q Q 



1-1 



fl 2 UJ 2 



n n 



dd log u 



2tt 2 1 



-Tr ( QoQ'' 1 
uj V 



(3.25) 



N ^ Tr ^ fi /_1 ^W _1 - J 1 (sW -1 ) T Jj J ■ (3.26) 

The boundary conditions for these master equations coincide with the weak coupling solutions 
(J33U). 



3.1.3 Energy density and flux densities 



Rewriting the energy density (13.91) in terms of our new variables and remembering the boundary 
term which vanishes when integrating over the entire plane, while nevertheless produces a big 
difference between the magnetic flux density and the energy density, we have 



£ = 2£dd\ogu + 2dd -TrfW 



r- 



However, the total energy 



E 



£ = 2tt£v 



2n£k 
n 



(3.27) 



(3.28) 



is simply proportional to the topological charge as always. 

The Abelian magnetic flux density is the first term (up to a factor) in the energy density 

(3.29) 



B = F® 2 = -2V2N ddlogu 



whereas the non- Abelian flux is the matrix defined in Eq. (I3.18p . The Abelian electric field 
density reads 



while the non-Abelian electric field density is given by 

A-jr 

E^t a = F^t a = —diTr (HHh a ) t a . 
A 4 

This can be written for G' = SU (N) as 



E°t a = —di 
H 



- [S'^W^S' - ^Tr ( 
u \ N \ 



(3.30) 



(3.31) 



(3.32) 



while for G' = SO(N) or G' = USp(2M) it is 



E?t a = -di 
3.2 Solutions 



^-s'- 1 m fi ,_1 - J ] (sw _1 ) T J^j s' 



(3.33) 



In the Abelian Chern-Simons theory, there exists a rigorous existence proof of the solutions in 
Ref. [34]. To our knowledge this has not rigorously been proved in the theory at hand. In the 
case of the vortices in the Yang-Mills-Higgs theory, the "covariant holomorphic" condition on the 
Higgs fields DH = 0, which is solved by the moduli matrix formalism, does uniquely determine 
the full moduli space of vortices via the Hitchin-Kobayashi correspondence [3!)ll3"oT I3TI 13"%] . which 
however has only been proved on compact spaces. This means that the corresponding master 
equations do not induce further moduli. For the vortices with the U(N) gauge group, an index 
theorem has been given in Ref. [9] while for generic gauge groups (under certain conditions) an 
index theorem has been given in Ref. [T3] . The index computed gives the number of moduli and 
does indeed correspond to the number of moduli found in the moduli matrix. 

The first part of constructing a solution is to write down the moduli matrix. Here we simply 
follow the way paved by the paper [12] using holomorphic invariants of the gauge subgroup G' . 
This boils down to some constraints for the moduli matrix to obey. A few examples of interest 
here is the case of G' = SU (N) 

det H {z) = z k + (z*- 1 ) , (3.34) 
while in the case of G' = SO, USp, respectively, we have 

m 2fe / 2k i \ 

H^(z)JH (z) = z~oj + [z^~ l \ , (3.35) 

where k is the vortex number (recall that v = ^ is the U(l) winding) and no = 2 in case of 
SO{2M) and USp(2M) while n = 1 for SO(2M + 1), M being positive integers. For SU(N), 
however no = N. 

The rather complicated looking master equations found in the last section are assumed to 
have a unique solution for each moduli matrix H (z) (up to V equivalence, see Ref. p2]). That is 
the moduli matrices are redundant and have to be identified by the following V transformation 

H (z) ~ V(z,z)H (z) , S(z,z) ~ V(z,z)S(z, z) , V G G c . (3.36) 

Here we conjecture the existence and uniqueness of the solutions to the master equation for each 
moduli matrix (up to the V equivalence). To provide plausibility for this claim we shall continue 
in two directions. 
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First we consider the weak coupling limit k — > and fi — > 0, which seems like an odd limit 
to take, but having an advantage. Looking at the theory (12.11) it is immediately seen that the 
matter fields are forced to stay in the vacuum manifold corresponding to the strong coupling 
limit of the normal non-Abelian vortex (i.e. with only a Yang-Mills kinetic term). In turn, this 
gives us a unique solution which in fact is the same solution as found in the strong coupling limit 
of the non-Abelian vortex with only a Yang-Mills kinetic term. This solution, appropriate only 
for vortices of the semi-local type, are usually called lumps in the literature. 

The second direction we will take will simply be to find some solutions by numerical calcula- 
tions. 

Now the existence of the solutions to the master equations, as we argue, makes it possible to 
exploit a lot of results developed in the literature. In short, 

the moduli space of non-Abelian Chern- Simons k vortices with gauge group G 
is equal to the moduli space of the non-Abelian Yang- Mills k vortices with gauge group G. 

(3.37) 

Moduli spaces of the non-Abelian vortices in N = 2 sQCD has been found in the literature in 
Refs. 01ITII for U(N) and in Refs. [13] for SO(N),USp(2M). 

Here we will summarize a few results from the literature. In the pioneering papers [9j [TO] 
discovering the non-Abelian vortices with gauge group U(N) (in contrast to the formerly found 
'Ln strings) the moduli space of a single vortex string was found to be 

M k=1 , G ,=su(N) = C x CP N ~\ (3.38) 

where the first factor denotes the position in the transverse plane while the second factor are 
orientational modes. For well separated k vortices, the moduli space can be composed as simply 
the symmetric product of that of the single vortex. This is not the case, when the centers coincide. 
In the k = 2, U(2) case, the moduli space has been found explicitly in the Refs. [151 US] 

M k=2 , G >=u{2) = C x WCP$ tl>1 , (3.39) 

which decomposes into a center-of-mass position and a weighted complex projective space with 
unequal weights giving rise to a conical type of singularity. In Ref. [J3] the moduli spaces of 
vortices with gauge groups G = U(l) x SO(N) and G — U(l) x USp(2M) has been found. A 
complication arises due to the fact that already for Nf = N flavors, the vortices are in general of 
the semi-local type (i.e. they have polynomial tails in their profile functions). The spaces quoted 
here correspond to the vortices of local type, thus they are constrained to have holomorphic 
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invariants with coincident zeroes. In the language of Ref. [13] this is obtained by constraining 
the vortices by the so-called strong condition 

m 2fc 

H^(z)JH (z) = (z- z )^ J . (3.40) 
The single local vortex with G' = U Sp(2M) has the moduli space 

M k=1 , G '=us p{ 2M) = C x US J^ ) , (3.41) 
while in the case of G' = SO(2M) it is found to be 

^-■(«=«)«( c »W. ■ 



lk=l,G'=SO{2M) 

where the ± denotes the chirality as described in detail in Ref. [IB] which is deeply rooted in the 
fact that the first homotopy group has in addition to the integers a Z 2 factor. This can also be 
interpreted as two spinor representations which is exactly the irreducible representations of the 
dual group G', where the dual is defined as being the group having the root vectors a* = -^=. 
For the k — 2, G' — SO(2M) the following orientational moduli spaces have been found to be 
locally 

M k=2 , G/=SO (, m) , Q22 = + i = x x Z 2 , (3.43) 



ka _ irom-1 v SO (Am) 

M k=2 , G/=S o(, m) , Ql2 =-i - M + x ^ - USp{2)m _ 1 - so{2) , 



SO (Am + 2) 

.M fc=2iG , =so(4m+2)iQz2=+1 = R™ x u{i) ^ USp(2 J )m x Z 2 , (3.45) 



gQ(4m + 2) 
USp[2) m x SO (2) 



M fc=2iG , =so(4m+2)jQz -_! = M+ x rreT _ /o ; m • (3.46) 



In the case of k — 1, G' = SO(2M +1), the moduli spaces are quite similar to the k = 2 even 
case. 



3.3 Weak coupling limit 
3.3.1 G' = U(l) x SU(N) 

Taking k = \i — > 0, we obtain from the D term conditions 

1 iV 1 N 

n' = (det Q y* fio , w = -7- (det fi )* , ^ = -7-^0 , (3.47) 
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which can be packaged together as a U (N) field Q. Instead of taking both couplings simultane- 
ously to weak coupling, we can play a game of taking only one of them, keeping the other finite 
(non-infinitesimal). Taking k — > and keeping \x finite we obtain 

u = -TrQ fi /-1 , (3.48) 

at the zeroth order in k while at first order we get the constraint 

iVTr M^o^'" 1 ) J = (TrlW" 1 ) 2 . (3.49) 

We note that only the Abelian field is determined, however at first order in the coupling constant 
we obtain a single constraint on the non- Abelian fields. Taking instead \i — > keeping k finite 
we have 

Q' = AQ , with A e const. , (3.50) 
to both zeroth and first order in \x. 



3.3.2 G' = U(l) x SO(N) and G' = U(l) x USp(2M) 

Taking k = \i — > we have from the D term conditions [T3l [33] 

Q' = H (z) . N M{z) , u = IttVWM , (3.51) 

where M = Hq(z)JH (z) is the meson field of the SO, USp theories according to the choice of 
the gauge group and in turn invariant tensor. 

A comment in store is that the Chern-Simons term is simply switched off in this limit and the 
lumps are the same as the ones living in the Yang-Mills theories experiencing infinitely massive 
gauge bosons. The point here, however, is to argue by continuity the existence and uniqueness of 
the solutions to the master equations for a given moduli matrix H (z) (up to the ^-equivalence 
relation). 



3.4 Numerical solutions 
3.4.1 Example: U(N) 

Let us do a warm-up and consider the single U(N) Chern-Simons vortex (k = /i) as has been 
found in Refs. [251 El] , however doing it in our formalism. Taking a simple moduli matrix 

H (z) =diag(2,liv-i) , (3.52) 
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which of course satisfies the constraint (13. 34ft . thus we can use the Ansatz for Q 

Q = e i 'diag(e {N - 1)x ,e- x l N _ 1 ) , (3.53) 
leading to the two coupled equations of motion 

dd [J, + (N - 1) X ] = \z - zof e-^ N ~^ (\z - z \ 2 e-^ N ~^ - |) , (3.54) 

Bd [j, - x ] = (e-* + x - jj) . (3.55) 

Notice that the two equations decouple in the sense that there only appear the combinations 
ip + (N — l)x and ip ~X- m f ac t ^ is easily seen that in this case, the field combination ip — x can 
be in the vacuum in all C which trivially solves the second equation. However, the first equation 
still needs to be solved numerically. The boundary conditions are 

V^oo = log f ~^~ J ' = log (M^) ■ (3.56) 

The equations become essentially Abelian when the couplings are equal k — fi, as was noted in 
Ref. [25] • The energy density is given by 

£ = 2£ddijj + 2dd [\z\ 2 e-^- {N - 1)x + (N - l)e"^ +x ] , (3.57) 

where the last term is the boundary term which of course integrates to zero. The Abelian and 
non-Abelian magnetic flux densities are given by 

F° 2 = -2v / 2iV ddip , = -2y/2N(N-l) Bd X t , (3.58) 

where the following matrix has been defined for convenience 

" vW^f ^" 1 ' -1 "- 1 '' (3 ' 59) 

which is traceless and has the trace of its square normalized to one half. The Abelian electric 
field density is 



Er = —\[^d r r r 2 e -^-(A r -Dx + (iV _ l) e -<Hxl 
K V N J 



(3.60) 



while the non-Abelian electric field density is 

2?r 



EX = J2N{N - l)d r [ r 2 e -V'-(JV-i) x _ e -v+xl t _ ( 3-61 ) 

We will find in the next subsection, that the numerical solution for this vortex for N = 2 is up 
to rescaling of some parameters equivalent to the vortex studied in the next subsection (when 
k — fi). Thus the concrete graphs are shown only for the vortex solution below. 
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3.4.2 Example: U(l) x S0{2M) and U(l) x USp(2M) 
Let us take a simple example of a moduli matrix 

H (z) = diag (zl M , 1m) , (3.62) 
which surely satisfies the constraint (13.351) . We take the Ansatz 

fl' = diag (e x l M , e" x l M ) , w = , (3.63) 
where det f2' = 1 is manifest. The equations of motion in terms of the new fields are 
dd X = - ^ (\z\ 2 e-^x + - j^j (\z\ 2 e-^-* - e~^*) 



- ^ ((k|V^) 2 - (e-^) 2 ) , (3.64) 



- — (blV^ - e-^) 2 . (3.65) 

It is interesting to note that under rescaling of the FI parameter £ — > M£, the above equations of 
motion are exactly the ones of the U(l) x SU(2) theory with the Ansatz used in the last section. 
The boundary conditions are 

^oo = log f~|~l*U ' Xoo = log(l^l) , (3.66) 

and the energy density reads 

£ = 2£ddiP + 2Mdd[\z\ 2 e-^- x + e-^ +x ] , (3.67) 

where the first term is proportional to the Abelian magnetic flux density 

Ff 2 = -4VM ddi; , (3.68) 

and the last is the boundary term which integrates to zero, while the non-Abelian magnetic field 
density reads 

F*t a = if 2 A t = -4VM Bd X t , t = ^^diag (1 M , -1 m ) . (3.69) 
The Abelian electric field density reads 



£r = ^ r [ r V^ + e-^] , (3.70) 
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0.5 1 1.5 2 2.5 3 3.5 0.5 1 1.5 2 2.5 

r r 

(a) (b) 

Fig. 1: (a) Profile functions for three different values of the coupling constants; a: k = 4, fi = 2; 
b: k = 2,fi = 2; c: K — 1, ^ — 2; the functions are plotted in traditional style with the winding 
field rising linearly and the non-winding field being constant at the origin. The FI parameter 
£ = 2. Notice that the VEV for these functions is 2~5. (b) The energy density £ for the vortex 
for the same three different values of the couplings. All the energy densities integrate to 7r£, 
within an accuracy better than ~ 10 -4 . 



whereas the non-Abelian electric field density is 

E«t a = E? A t = ^^d r \r 2 e^-x - e^+*l t . (3.71) 
/i 

We show the vortex with this Ansatz corresponding to different values of the coupling constants 
ft, /i in the following figures. Here we will take for definiteness the group G' to be 5*0(4) or 
USp(4:) hence M = 2, which within the chosen Ansatz are equivalent. We furthermore set £ = 2. 
The total energy is thus (recall the Ansatz is for a single k = 1 vortex) 

E = [ £ = tt£ . (3.72) 
Jc 

In Fig. we show the profile functions of the vortex in the traditional way, where the color-flavor 
matrix is parametrized as follows 

= diag(/(r)e ie l 2 ,s(r)l 2 ) , (3.73) 

which of course is equivalent to the parametrization in terms of ip,X- in Fig- [D 3 t ne energy 
density of Eq. (13.671) is shown. The integral of the energy density is identically equal to the 
integral of the Abelian magnetic flux, as it should be. We see the vortex size is proportional to 
the coupling constants. In Fig. [2] we show the Abelian (a) and the non-Abelian (b) magnetic 
field, respectively. We observe that the Abelian magnetic field is negative at the origin while 
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0.5 1 1.5 2 2.5 3 3.5 4 0.5 1 1.5 2 2.5 3 

r r 

(a) (b) 

Fig. 2: (a) The Abelian magnetic field F® 2 (trace-part) for three different values of the couplings. 
Notice the equal coupling case has zero magnetic field at the origin while the different coupling 
cases have negative and positive values, respectively, (b) The non-Abelian magnetic field F^ 2 
(traceless part) for different values of the couplings. Notice the opposite behavior of the non- 
Abelian magnetic field with respect the Abelian one at the origin, see also Fig. |H The FI 
parameter £ = 2. 



the non-Abelian magnetic field is positive, in the k = 4, /i = 2 case. The contrary holds in the 
k — 1, fi = 2 case where the non-Abelian magnetic field is negative at the origin while the Abelian 
field is positive. It turns out that the combination 



(kF° 2 + / < 2 a )| =o. 



(3.74) 



An immediate consequence is that for \k\ ^> \F^ 2 A \ ^> \F® 2 \ at the origin and vice versa. Plots 



of the Abelian and non-Abelian magnetic fields normalized as in Eq. (13.741) are shown in Fig. H] 
with k — 4, fi — 2 in (a) and k = = 2 in (b), respectively. At the origin this combination 
cancels to a numerical accuracy better than 10 -5 . First let us demonstrate the formula (13.741) by 
calculating the fields in the limit r — > 



kf: 



I 

12| r ->0 



-^1 N 2 A | 



TWO 



AttVM 



K 



—)+- (e"' 
M) fx y 



(3.75) 



Note that the value of the magnetic fields only depends on the field combination ip — x, and 
it is understood that it has to be evaluated at the origin in the above equation. Secondly, let 
us demonstrate that the magnetic fields are zero at the origin in the case of equal couplings. 



15 



0.5 1 1.5 2 2.5 3 3.5 0.5 1 1.5 2 2.5 3 

r r 

(a) (b) 

Fig. 3: (a) The Abelian electric field in the radial direction E r (trace-part) for three different 
values of the couplings, (b) The non- Abelian electric field in the radial direction E^ A (traceless 
part). The FI parameter £ = 2. 



Subtracting Eq. (13 .64ft from Eq. (13. 65ft we have 



dd (ij> - x) 



7T 



K 



K 



\z\ 2 e^-x 



) 2 + i 



2e 



M 



\z\ 2 e-*-* 



+ 1 + 



■>P+x\ 2 _ 



M 



-ll>+X 



(3.76) 



which depends on z, z when the coupling constants are different, k ^ \x. However, when the 
coupling constants are equal, Eq. (I3.76P reads 

4tt 2 ( i 



dd (ip - x) 



-4>+x 



2M 



(3.77) 



which allows the field combination ip — x to stay constant with the value 

;iM y 

i'-x = log 



(3.78) 



Plugging this (constant) solution into Eq. (13.751) we obtain readily F^ 2 = F^ A = in the limit 



0. 



In Fig. [3] is shown the Abelian (a) and non- Abelian (b) electric fields with different values of 
the couplings. 

In Fig. [5] we show a sketch of the magnetic fields of Abelian and non- Abelian kinds, respec- 
tively, in the case of k > fi (a) and in the case of k < /i (b). The integral over the plane of the 
Abelian magnetic field density is proportional to the topological charge of the vortex, the winding 
number which in turn renders the soliton topologically stable. The vortex solution with negative 
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K =4,n=2 k=1,h=2 




(a) (b) 

Fig. 4: Differently normalized Abelian and non-Abelian magnetic fields as kF® 2 and fiF^ 2 A for 
(a) k = 4, /i = 2 and (b) n = = 2. This combination cancels exactly at the origin (to a 
numerical accuracy better than ~ 10 -5 ). The FI parameter £ = 2. 




Abelian magnetic field Non-Abelian magnetic field Abelian magnetic field Non-Abelian magnetic field 



(a) k > jj, (b) k < /x 

Fig. 5: (a) Sketch of the magnetic fields where the Abelian (red/left) is negative at the origin 
and the non-Abelian (blue/right) is positive for k > \i. (b) Contrarily the Abelian (red/left) is 
positive at the origin while the non-Abelian (blue/right) is negative for k < /i. 
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0.5 1 1.5 2 2.5 3 3.5 0.2 0.4 0.6 0.8 1 1.2 1.4 

r r 

(a) (b) 

Fig. 6: (a) Profile functions for three different values of the coupling constants; a: k = —4, fi = 2; 
b: k = — 2,/i = 2; c: k = — = 2; the functions are plotted in traditional style with the 
winding field rising linearly and the non-winding field being constant at the origin. The FI 
parameter £ = 2. Notice that the VEV for these functions is 2~5. (b) The energy density £ for 
the vortex for the same three different values of the couplings with opposite signs. All the energy 
densities integrate to 7r£, within an accuracy better than ~ 10" 4 . Notice that the extrema of 
the energy density is at the origin, just as in the case of the ANO vortices or the non-Abelian 
generalizations. 



winding number k < can be interpreted as an anti-vortex. Hence, one could wonder which 
interpretation to give the small substructure found in this vortex solution - a small anti-vortex 
trapped in the non-Abelian vortex, as a bound state, not rendering the solution unstable. 

Opposite signs of coupling constants 

We will now consider taking one of the couplings to be negative, say n < and fi > 0. Choosing 
both signs negative yields the same solution as already mentioned, however with flipped electric 
fields. In the case of k > and fi < 0, the solutions are equivalent to the ones we will consider 
now, just with the signs flipped of the electric fields. The Chern-Simons characteristics have 
been lost in this case, the vortex instead has the magnetic field concentrated at the origin - just 
as in the case of the ANO vortex or the single U(N) non-Abelian generalization. In Fig. [HI the 
profile functions and energy densities for different solution are shown. In Fig.[7]the corresponding 
magnetic fields are shown while in Fig. [8] the electric fields are shown. 
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Fig. 7: (a) The Abelian magnetic field Ff 2 (trace-part) and (b) the non-Abelian magnetic field 
F^ 2 A (traceless part) for three different values of the couplings with opposite signs. Notice that 
the magnetic field density resembles that of the ANO vortex or the non-Abelian generalizations, 
viz. they have the extrema at the origin. The FI parameter £ = 2. 




0.5 1 1.5 2 2.5 3 3.5 0.5 1 1.5 2 2.5 



r r 

(a) (b) 

Fig. 8: (a) The Abelian electric field in the radial direction E r (trace-part) and (b) The non- 
Abelian electric field E^ A (traceless part) for three different values of the couplings with opposite 
signs. The FI parameter £ = 2. Note that the electric fields are back-to-back. 
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4 Discussion 



We have thus brought the powerful moduli matrix formalism into the non-Abelian Chern-Simons 
model (which supports topological non-Abelian vortices), and have conjectured that the moduli 
spaces of the non-Abelian vortex solutions of these systems are indeed identical to those of the 
vortex solutions in the Yang-Mills-Higgs models with corresponding gauge groups. We have not 
proved that every moduli matrix has a unique and existing solution to the master equations found. 
Nevertheless we have argued the plausibility of such a claim by taking the weak coupling limit 
which immediately yields the lumps of the Yang-Mills-Higgs models, as it is just the algebraic 
solutions to the D term conditions. 

We have then studied some numerical solutions of non-Abelian vortices, by choosing an Ansatz 
to the master equations, working mainly with the G' = 50(4) and G' = USp(4) gauge groups. 
We have studied the case of different couplings with both couplings positive yielding vortex 
solutions with a small negative Abelian (non-Abelian) magnetic field density at the origin and a 
corresponding positive non-Abelian (Abelian) magnetic field density, which have a combination 
that is always zero (at the origin). Keeping the couplings equal provides the typical Chern-Simons 
characteristic that the magnetic field vanishes at the origin yielding a ring structure. These new 
type of solutions could perhaps be interpreted as an anti-vortex sitting inside the non-Abelian 
vortex as a stable bound state, with the stability provided by topological arguments. 

An interesting question is to which extent this substructure found in the non-Abelian vortex 
solutions alters the dynamics of the vortices. 

Furthermore, by changing the relative sign of the coupling constants a vortex solution with 
the magnetic field density concentrated at the origin has been found. 

An obvious future study related to these vortices and also to the ones of Ref. [13] could be 
to make an explicit construction with exceptional groups and investigating the corresponding 
moduli spaces. Especially interesting would be the center-less groups. 

Another interesting path to follow is to consider the construction of the non-Abelian vortices 
in Chern-Simons models with more supersymmetries, e.g. considering the model of Aharony- 
Bergman- Jafferis-Maldacena [39] . An Abelian non-relativistic Jackiw-Pi vortex has already been 
found in this model [40J. Another attempt to construct vortices in the latter model has recently 
been made, resulting in the non-Abelian vortex equations of the Yang-Mills-Higgs models |41j . 
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